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Abstract 

Consider the random matrix obtained from the adjacency matrix 
of a random <i-regular graph by multiplying every entry by a random 
sign. The largest eigenvalue converges, after proper scaling, to the 
Tracy-Widom distribution. 

1 The result 

Let G = (V,E) be an element of $(N, d), i.e. a random uniformly chosen 
(i-regular graph with N vertices, and let S : E — > { — 1, +1} be independent 
random signs: 

¥{S(u,v) = 1} = ¥{S(u,v) = -1} = 1/2 . 
Consider the N x N (random) Hermitian matrix M, 




(u,v) e E 
(u,v)<£E 



with eigenvalues Ai < A2 < • • ■ < \n- 
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Theorem. Assume that M is as above, with 3 < d <C N 2 ^ 3 . The scaled 
largest eigenvalue 



converges in distribution to the Tracy-Widom law TW\. The same is true 
for the scaled smallest eigenvalue 



The proof is based on the combinatorial arguments from [2], which can be 
seen as a modification of Soshnikov's approach [9] to the asymptotic distri- 
bution of the extreme eigenvalues of random Hermitian matrices with inde- 
pendent entries. Another crucial ingredient is an estimate of McKay [I] on 
the number of subgraphs of a random graph. 

2 Connection to non-backtracking walks 

Consider the matrices 




(1) 





M (n) = (d- l) n/2 { U, 



n 



M 



d-1 



1 



C/n-2 




where 




are the Chebyshev polynomials of the second kind, and formally 



17-2 = U. 



. 



It is not hard to see (cf. [S]) that, for any u ,u n G V, 




where the sum is over all paths p n = u uiu 2 



u n such that 



(a) (uj, Uj + i) G E for any < j < n — 1; 
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(b) Uj +2 7^ Uj for < j < n — 2 (p n is non-backtracking.) 
Therefore 

tr Af< n) = ^"M ttottl M UlM2 • • • M Un „ lU „ , 
where the sum is over paths p n satisfying (a), (b), and 

(c) u n = Mo- 
Taking the expectation (over S) and noting that EgM^, = ((— l) k + l)/2, we 
deduce: 

Lemma 1. TTie expectation E^-trMW zs egna/ to i/ie number P n (G) of paths 
p n that satisfy (a), (b), (c), and 

(d) even/ edge (u, i>) G -E appears an even number of times on p n . 

In particular, E tr M^ n+l ^ = 0, we shall therefore study E tr = E G P 2n (G) , 

where E^ denotes expectation over the random choice of the graph G. 

3 Diagrams 

Following [2], we group the paths p 2n satisfying (a)-(d) into (topological) 
equivalence classes, which are in one-to-one correspondence with diagrams: 

Definition 2. Let (3 G {1,2}. 

• A diagram is an (undirected) multigraph G = (V,E), together with a 
circuit p = u Ui ■ ■ ■ u on G, such that 

— p is non-backtracking (meaning that no edge is followed by its 
reverse, unless the edge is a loop uu); 

— For every (u, v) G E, 

# {j I uj = u, u j+1 =v} + #{j\ uj = v, u j+ i = u} = 2 ; 

— the degree of Uq in G is 1; the degrees of all the other vertices are 
equal to 3. 

• A weighted diagram is a diagram G together with a weight function 
w {-1,0,1,2,- ■■}. 
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Informally, if (ui, u 2 ), (u 2 , u 3 ), • • • , (iik-i, Uk) is a sequence of edges with 

degwi = degw fc = 3, degw 2 = • • • = degw fc _! = 2 , 

we replace U\U 2 • • • with a single edge e = u k ), to which we assign the 
weight w(e) = k — 2. Thus, the path 



corresponds to the diagram in Figure 1 (left), whereas 

1.2.3.4.5.6.7.8.9.7.8.9.6.5.10.11.12.13.14.15.16.12.13.14.15.16.12.11.10.5.4.3.2.1 
corresponds to Figure 1 (right). 



Figure 1: Some diagrams: s — 1 (left), s = 2 (center, right) 

More complicated paths (with vertices of degree greater than 3) also cor- 
respond to weighted diagrams; the formal correspondence is described in [21 



The following lemma summarises the results of [21 II.1TI.2] that we shall use: 
Lemma 3. 

1. (Claim II. 2. 2.) All the diagrams can be classified according to a pa- 
rameter s = 1,2, ■■ • (see Figure^) A diagram with parameter s has 
#V = 2s vertices and #E = 3s — 1 edges. 

2. (Proposition II. 2. 3) The number D\(s) of diagrams corresponding to a 
given value s satisfies 



1.2.3.4.5.6.7.4.6.5.7.4.3.2.1 



o 




III]- 



(s/Cf < Di(s) < c s -\s s , 



where C > 1 is a numerical constant. 
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3. For a diagram 5 with parameter s, let W(S) be the number of associ- 
ated weighted diagrams, and W(8) - the number of associated weighted 
diagrams with positive weights. Then 



;;!;)<-m<-w {s) <(^_- 2 



In particular, if s = o(y / n), 

W{5),W(6) = (l + o(l)) 



n 3s 2 



(3s - 2) 



4- Every path corresponding to 5 has at most n — s + 1 distinct vertices. 
If the associated weighted diagram has positive weights, the path has 
exactly n — s + 1 distinct vertices; of these, 1 is of degree 1, n — 3s + 1 
are of degree 2, and 2s — 1 are of degree 3. A general path is obtained 
from one as above by identifying several pairs of connected vertices, and 
erasing the connecting edges. 

5. (Claim II. 1. 4) The number of different paths corresponding to a given 
weighted diagram is at most jV n ~ s+1 . // the weights are positive, the 
number is exactly 

N(N - i)...(N-n + s) . 

4 Counting the paths 

Fix a diagram 5 with parameter s. We shall now estimate the expected 
number of paths p 2n associated to 5 in a random graph G. We shall use the 
following special case of jH Theorem 2.10]. 

Let L be a subgraph of the complete (labelled) graph K N ; let £ u be the 
degree of u in L, E G = #E(G) = dN/2, and E L = #E(L) = £4/2- 
Finally, denote = a!/(a — b)\ 

Proposition 4 (McKay). If E G — E L > 3d(d + 1), then the probability F L 
that a random d-regular graph G contains L satisfies 



n n d [iu] 
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where 

< * d(d+l) ) El [e 

2(E G -E L -2d(d+l)) I 



i + f . [ (Eq-iw 



(E G -2d 2 Y E ^ ' 
We shall deduce the following estimates: 
Lemma 5. 

1. Forn < N, 



m n (G) <n(d-2)(d-l) n - 1 expl 1 ^ w (l + d/V^N)\ , (4) 



where C > is a numerical constant. 
2. For d 2 /N < n < min( v /n, N/d), 



, " \ - 

i d s ~ l (d-l) s N s - 1 (3s -2)! 



(5) 



Proof. Let p2 n be a path associated to a diagram 5 with parameter s, and 
let L be the induced (undirected) graph: 

E{L) = {{ Uj ,u j+1 )} . 
Then, by item H] of Lemma [3j 

Yl u d [£u] d n - s+1 {d - l) n ~ s {d - 2) 2s - 1 



2 e lE I e l] ~ 2 n (d/V/2)M 

G 



{d-2) 2s ~ l (d-l\ n {dN/2) n ^ 



~ d s ~ l (d-l) s \ N J (d/V/2)M 

(CW\ fd-V\ n (d-2) 2s - 1 
- eXP \ dN J V N ) d-i(d-l) 
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where C% > is a numerical constant. Also, 

E(L) = (Q?iY/2)[ "' < exp {C 2 nd/iV} ; 
V ; (dN/2-2d 2 )^ ~ P1 1 5 

therefore 

Cm 2 C 2 nd\ (d-l\ n (d-2) 23 - 1 



Fl < exp 



dN N J V N J d'-^d-l 



According to item [3J the number of ways to choose the weights on 5 is at 
most 

'n + 3s - 2\ n 3s - 2 
3s - 2 J - (3s -2)! ' 

and, according to item HJ the number of ways to choose the vertices is at 
most jV n ~ s+1 . Finally, the number of diagrams with parameter s is at most 
C s ~ 1 s s . Therefore 



rn-s+l n 



(3,-2)! 

fCX C 2 ndl /d-l\ n (rf-2) 2s - 

x exp < — h 



dN N } \ N J d'-^d- l) s 
< n(d -2)(d -1)" exp {^ + ^} (7) 



S>1 



(2(^-1))! 

< n ( d _ 2 ){d - l)"- 1 exp + d/V^N) 

This proves the first statement. 

Suppose d 2 /N <Cn< mm(\/~N, N/d). Choose so such that 

's = l, n<N 1 ^ 
n a/2/Ni/2 < So < ^ , n > N l ' A , 

For s < s and positive weights, the inequalities above are asymptotic equal- 
ities. Namely, 

'd-l\ n (d-2) 2 *- 1 



:i + 0(1)) > N J d-Kd-1) 
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W(8) = (l + o(l))- 



n 3s 2 



(3s-2)! ' 

Hence the contribution of paths with positive weights and s < so is 



1<S<S| 

The contribution of paths with s < Sq and not all weights positive is negligi- 
ble, since 

W(S) - W(S) = o(l)^^ . 

The contribution of paths with s > sq is also negligible, by the estimates in 
the proof of item 1 of this lemma (see (J7|) above.) Thus 



EP 2n (G) = (l + o(l)) D ^ Nn ~ S+1 ji^V. ( 

1<S<S ^ ^ 



rf-lV (d-2] 



2s-l 



N J d s - l (d-l) s 

Ks<oo 



□ 



5 Conclusion of the proof 

Applying Lemma[T]and using the definition of M^ 2n \ we immediately deduce: 
Corollary 6. 

1. Forn < N, 

Etrf/ 2n j^-|L=} <nexp|^^(l + d/v / ^V)| , 

where C > is a numerical constant. 

2. For d 2 /N < n < min( v / iV, iV/d), 

_ rr f Ml , . , . ^ /(d-2) 2 n 3 Y _1 D a (s) 

Etr M^H 1+o(1))n i;(Mw) ps^ji- 
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Now, by [21 Theorem 1.5.3, Proposition II.2.5], an N'xN' random matrix A in 
from the Gaussian Orthogonal Ensemble GOE(iV') satisfies, for n = o(viV') 



EC/ 2n (A inv /(2v / iV 7 ^2)) = (1 + o(l)) n ^(n 3 /iV) 

S>1 



DM 
(3s - 2) 



Hence, for N' = [ ^ 2 ^ A^j , the second item of Corollary [H] implies: 



EtrU- 



2n 



j M 

\2y/J=T 



(l + o(l))EtrC/ : 



2» 



J ^inv 



Applying the arguments of [21 II. 3], this is extended to expectations of prod- 
ucts of several such traces. 

As in [21 1.5] this implies that the random point process 



d(d 



(d-2y 



N 



2/3 r 



A, 



2Vd 



N 



converges (in the appropriate topology, [21 1.2]) to the Airy point process 2Ui, 
and hence the distribution of ([!]) converges to TW\. 

□ 



6 Some remarks 

1. This note continues the study of the distribution of the scaled largest 
eigenvalues of random matrices with independent entries. For the Gaussian 
Orthogonal Ensemble, the limiting distribution was identified by Tracy and 
Widom, building on several earlier works; see [10] and references therein. 
This result was extended to a rather general class of matrices with indepen- 
dent entries by Soshnikov [9J. Important further results have been obtained 
by Ruzmaikina [7] , and Khorunzhiy and Vengerovsky [3] . A different proof of 
Soshnikov's theorem was given in [2]; it is the basis of the present argument. 

2. In the physical literature, the spectrum of random matrices similar to 
M has been recently studied by Oren, Godel, and Smilansky [6], who have 
focused on the bulk of the spectrum (rather than the extreme eigenvalues). 

3. The conclusion of the theorem remains valid for more general weights 
S(u,v), such that 
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(Al) The distribution of every S(u,v) is symmetric; 
(A2) ES(u,v) 2k < (C k) k ; 
(A3i) ES(u,v) 2 = 1. 

The extension of the proof follows the lines of [21 Part III]. 

4. One may also consider complex weights S(u,v) that satisfy (Al), (A2), 
and 

(A3 2 ) E($S{u,v)) 2 = E{%S{u,v)) 2 = 1/2, E{^S{u,v)QS{u,v)) = 

(for example, S(u, v) may be uniformly distributed on the unit circle.) Then 
the scaled extreme eigenvalues ([]]), (j2J) converge in distribution to the Tracy- 
Widom law TW 2 . 

5. If all the weights S(u,v) are equal to 1, the matrix M is the adjacency 
matrix A G of G. Obviously, the largest eigenvalue of A G is equal to d. 
Numerical evidence due to Miller, Novikoff, and Sabelli [5J indicates that the 
second largest eigenvalue A 2 of A G converges to TW\ after proper scaling. If 
true, this would have important consequences, see [Sj and references therein. 

6. Our result can be reformulated as a bound for the "new" eigenvalues of a 
random 2-lift of a random c?-regular graph. We refer the reader to the work 
of Bilu and Linial pQ for the definitions. 

Acknowledgment. I am grateful to my supervisor Vitali Milman for his 
support and encouragement. I thank Nathan Linial, Steven J. Miller, Idan 
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